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Abstract 

We report electrical measurements of the current-induced spin polarization of the surface current 
in topological insulator devices where contributions from bulk and surface conduction can be disen¬ 
tangled by electrical gating. The devices use a ferromagnetic tunnel junction (permalloy/ALCH) as 
a spin detector on a back-gated (Bi,Sb) 2 Te 3 channel. We observe hysteretic voltage signals as the 
magnetization of the detector ferromagnet is switched parallel or anti-parallel to the spin polariza¬ 
tion of the surface current. The amplitude of the detected voltage change is linearly proportional 
to the applied DC bias current in the (Bi,Sb) 2 Te 3 channel. As the chemical potential is tuned 
from the bulk bands into the surface state band, we observe an enhancement of the spin-dependent 
voltages up to 300% within the range of the electrostatic gating. Using a simple model, we extract 
the spin polarization near charge neutrality (i.e. the Dirac point). 
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I. INTRODUCTION 


The combination of strong spin-orbit coupling, time-reversal symmetry and inversion 
symmetry is known to create gapless helical, Dirac surface states that lie within the bulk 
band gap of narrow band gap semiconductors such as the Bi-chalcogenides. These three- 
dimensional ( 3 D) “topological insulators ” 1 3 have begun to attract significant contemporary 
interest for spintronics since the spin-momentum locking in the surface state endows an 
inherent spin polarization to surface state charge currents ! 4 5 This has motivated several 
recent experiments on “topological spintronics.” Spin-transfer torque has been demonstrated 
in bilayers of ferromagnets and topological insulators,®^ with record values of the spin torque 
ratio at room temperature in Bi 2 Se3%nd at low temperature in (Bi,Sb) 2 Te3 . 7 Spin Hall angles 
have also been measured through the inverse spin Hall effect in topological insulators, created 
by the injection of a spin current from a metallic ferromagnet through spin pumping 9 ^nd 
through spin-polarized tunneling . 12 

The spin polarization of the surface states of topological insulators has been measured in 
extensive spin- and angle-resolved photoemission spectroscopy studies ^ 13 16 However, direct 
electrical measurements, relevant for developing spintronic devices, are just beginning to 
emerge and derive from a well-established methodology in spintronics: the essential idea 
is to use a ferromagnetic contact as a voltage probe of the spin-dependent electrochemi¬ 
cal potential. As the relative orientation of the current spin polarization and ferromag¬ 
net magnetization changes, so does the measured voltage. 1 ' This scheme was first used to 
clearly demonstrate spin-momentum locking using ferromagnetic tunnel contacts on Bi 2 Se3 
transport channels, albeit with the chemical potential of Bi 2 Se3 in the bulk conduction 
band and at temperatures below 150 K . 18 More recent experiments have extended such 
measurements to room temperature in Bi 2 Se3 19 and to topological insulators with reduced 
bulk conduction . 20 22 Further, a recent study used spin-polarized tunneling in ferromag¬ 
net/insulator/topological insulator junctions to map out the energy-dependence of the ef¬ 
fective spin polarization in BioSesi 12 Here, we demonstrate an experiment that is comple¬ 
mentary to these prior studies. We show that the standard three probe spin potentiometry 
scheme using a magnetic tunnel junction (MTJ) 171 ® can be effectively applied to electri¬ 
cally gated topological insulator devices, thus yielding insights into the chemical potential 
dependence of the spin polarization of the surface state current. In our measurements, we 
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map out the spin-dependent voltage as the chemical potential is tuned from the bulk con¬ 
duction band through the Dirac point and into the valence bulk band. We find that the 
spin-dependent voltages are enhanced as the chemical potential is moved away from the bulk 
conduct ion/valence bands into the surface state band. We also show that we do not observe 
possible contributions of the opposite spin polarization from Rashba spin-split, bands within 
the range of the electrostatic gating. 


II. SAMPLE SYNTHESIS AND DEVICE FABRICATION 

(Bi,Sb) 2 Te 3 thin films were grown on SrTi 03 substrates by molecular beam epitaxy 
(MBE ). 23 We used a Bi:Sb ratio of 1 to place the chemical potential into the bulk 
band gap. The samples had film thicknesses of 7 nm (SI), 8 nm (S2), and 6 nm (S3); this 
minimizes contributions from unpolarized bulk carrier conduction and allows the tuning of 
the electron chemical potential of both the bottom and top surfaces by back gating. We note 
that the film thicknesses are above the limit where Dirac point is disrupted by hybridization 
of the top and bottom surface si at e.Jh ' 21 251. 

For clean measurements of the surface state spin polarization, we developed a device 
geometry that restricts all measurements to the top surface [Fig|l](a)]. We first define a cross¬ 
shaped (Bi,Sb) 2 Te 3 channel by standard photolithography followed by Ar plasma etching. 
A bilayer photoresist window of circular or rectangular shape is then patterned at the center 
of the cross-shaped channel. The bare surface of the topological insulator is cleaned using a 
low-power Ar plasma etch (15 W for 20 s) that removes possible native oxides and photoresist 
residues along with ~1 nm (Bi,Sb) 2 Te 3 . Following the sputter deposition of an AI 2 O 3 seed 
layer (~0.3 nm), we use atomic layer deposition (ALD) for conformal deposition of an AFO 3 
layer ( 2.1 nm for devices on SI and 1.8 nm for devices on S 2 ) over the topological insulator 
surface. A Py layer (20 nm) and a thin Au capping layer (5 nm) are then deposited by e- 
bearn evaporation, followed by a lift-off process. Finally, the top of the MTJ is extended to 
two Au contact pads over a 60-nm-thick Si 3 N 4 layer to isolate the MTJ from the topological 
insulator channel as shown in Figs jljb) and[l](c). Each sample has multiple devices with 
various MTJ areas. The five devices studied in this paper are listed in Table I. 

For microstructural analysis, a device was cross-sectioned using a focused ion beam (FIB, 
FEI Quanta 200 3D) and then imaged using high-angle annular dark-held scanning trans- 
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mission electron microscopy (HAADF-STEM) [Figs jljd) and[l](e)]. For STEM imaging, an 
aberration-corrected (CEOS DCOR probe corrector) FEI Titan G2 60-300 S/TEM equipped 
with a Schottky X-FEG gun was operated at 200 kV with a probe convergence angle of 16 
mrad. Additionally, the AI 2 O 3 layer was identified as being amorphous by conventional 
TEM imaging (not shown), using an FEI Tecnai G2 F30 at 300 kV. 

III. RESULTS AND DISCUSSIONS 

A. Ambipolar transport of (Bi,Sb)oTe3 

We first discuss the back-gate-voltage (Vg) dependence of the channel resistance and the 
Hall resistance. The former is measured through the two-terminal channel voltage between 
leads 2 and 4 and the latter between 3 and 5 is measured while a DC current of 1 /iA flows 
between leads 2 and 4 [see Fig§ a) for lead configuration]. The channel between two Ti/Au 
leads (between 2 and 4 as well as between 3 and 5) is 300 pm long by 200 pm wide, and 
the contact resistance of each lead is ~100 D. The Hall voltage was measured as a function 
of applied perpendicular-to-plane magnetic field of 10 kOe. The results show the typical 
signatures of ambipolar transport in a topological insulator: the change of carrier type from 
electrons to holes as well as a peak in channel resistance near the charge neutrality point (at 
Vg — 20 V for SI and at Vg = 10 V for S2) [Figs|2](a) and[2](b)]. Carrier densities at each 
end of the applied gate-voltage range were obtained to be n e = 3.53 x 10 13 (2.44 x 10 13 ) 
cm -2 at Vg = 140 V and n p = 1.91 x 10 14 (1.69 x 10 14 ) cm -2 at Vg = —140 V for SI 
(S2), respectively [Fig|2](c)]. These values give an idea of the range of the chemical potential 
tuning. We believe that we were able to change the chemical potential across the bulk band 
gap from a position near the bottom of the conduction band (at 140 V) to a position below 
the top of the valence band (at —140 V). 

B. Weak anti-localization of (Bi,Sb)2Te3 

To further understand the surface and bulk conduction in the (Bi,Sb) 2 Te 3 films with 
respect to the gate voltage, we carried out magneto-transport measurements using the same 
wiring configuration. Near 2 K, the magneto-resistance of a topological insulator channel 
shows a sharp cusp near zero magnetic field as shown in the inset to the Figj2](d). This 
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suppression of resistance near zero magnetic field can be explained by weak anti-localization, 


typically seen in conventional 2D metals with strong spin-orbit coupling 26 ' as well as in 3D 
topological insulators that have a non-trivial n Berry’s phase 13 . The quantum corrections 
to the diffusive magneto-conductance are given by 27 



( 1 ) 


Here, a is the 2D channel conductivity and can be determined from the 2D channel resistivity 
p xx and 2D Hall resistivity p xy via a = p xx /(p xx + p xy ). Also, e, h, a, and 1$ are the 

electronic charge, the Planck’s constant, the digamma function, the pre-factor, and the 
coherence length, respectively. Different types of localization behaviors can be identified 
by the values of the pre-factor a, which takes on values of 1, 0, —1/2 for the orthogonal, 
unitary and symplectic classes, respectively] 2 ' The weak anti-localization behavior of the 
surface state of a 3D topological insulator is in the symplectic class with a = — 1/2 expected 
for one channel conduction of a surface state or bulk-carrier-coupled top and bottom surfaces. 

For two conduction channels of decoupled top and bottom surfaces, a = (—1/2)+ (—1/2) = 

_ 1 |28H30] 

Figure |2](d) shows the gate-voltage dependence of the pre-factor a obtained from the 
fitting of weak anti-localization by Eq.(l) at different gate voltages. At each end of the 
applied gate-voltage range (—140 V and 140 V), a approaches —1/2, indicating that the 
chemical potential is on or near the bulk bands whose carriers couple the top and bottom 
surfaces. Near the charge neutrality point where the bulk conduction is minimal, a deviates 
from —1/2 and approaches —1, indicating that the surface states on the top and bottom 
surfaces have become decoupled, independent conduction channels, a changes faster with 
the chemical potential as it is lowered from the charge neutrality point towards the valence 
band because the Dirac point of (Bi,Sb) 2 Te 3 is closer to the top of the valence band than 
the bottom of the conduction band. We note that the contact resistance was subtracted 
before fitting the two-terminal magneto-resistance data to Eq.(l). 

C. Tunnel junction properties 

The RA product of the tunnel junctions is in the range 10 6 -10' Mfh/irn 2 . We observe qual¬ 
itatively similar behavior of the gate-voltage dependence, the tunnel junction characteristics, 
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and the spin-dependent voltages from all the devices. Typically, the temperature dependence 
of the zero-bias junction resistance shows an insulating behavior [Figj3|a)], while the chan¬ 
nel resistance of the topological insulator shows metallic behavior below 100 K [Figj3^a) 
inset]. The junction resistance as a function of measured DC voltage ( V B c ) across an MTJ 
shows an anomaly at zero-bias [Figj3^b)] , commonly seen in MT.I.v* 1 . We also observed 
an asymmetry in the junction resistance between positive and negative voltages [Figj3]]b)]. 
Interestingly, the asymmetric behavior qualitatively follows the gate-voltage dependence of 
the channel resistance [Figj2j]a)] , which can be interpreted as probing of the density of states 
of the topological insulator top surface via tunneling between the Vnc-tuned metal and the 
Vd-tuned topological insulator. 

D. Spin-dependent voltages from topological insulator spin polarization 

Our three-terminal potentiometric measurement scheme follows the proposal by Hong et 
alW- We use an A^Oa/Py junction for the ferromagnetic electrode while a dc current flows 
through a topological insulator channel with an in-plane magnetic held perpendicular to the 
current direction. When the direction of the Py magnetization switches from parallel to anti¬ 
parallel to the current-induced spin polarization of the topological insulator surface state, 
we observe step-like changes in the detected voltage, resulting in hysteretic spin-dependent 
voltages associated with the coercive held (~100 Oe) of the Py layer. For positive currents 
as shown in Figsj4](a) and[4](b), the relative direction between the Py magnetization and 
the topological insulator spin polarization determines the positive change in the detected 
voltage AV = [V(M) — V{— M)]. For reversed currents [FigsQc) and[4|d)], the direction of 
topological insulator spin polarization is opposite to the case of positive currents, so that the 
change in the detected voltage AV becomes negative. Figures |4|e) and[4]]f) show that the 
measured voltage change is linearly proportional to the bias currents. This linear relation 
can be expressed by 

AV = IR B (p-m), (2) 

where R B , p, and m represent ballistic resistance, degree of the spin polarization per unit 
current in the topological insulator channel, and effective magnetic polarization of ferro- 
rnagnet, respectively. 1 ' Depending on the position of the chemical potential tuned by the 
electrical gating, the measured spin-dependent voltages may originate not only in the topo- 
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logical insulator surface states but also in the Rashba spin-split bands. Equation (2) is 
applicable to the spin-dependent voltages from combined conduction channels including the 
surface states, the Rashba spin-split bands, and the unpolarized bulk bands in both ballistic 
and diffusive limits. 


E. Gate-voltage dependence of spin-dependent voltages 

We now measure the hysteretic spin-dependent voltages while changing the chemical 
potential across the bulk bandgap. When the chemical potential is placed in the bulk 
band gap, the observed spin-dependent voltages are two to three times larger. Using two 
devices (Dev3 from SI and Dev4 from S2), multiple measurements of magnetic held sweeps, 
repeated at least 8 times, were carried out every 10 V from 140 V to —140 V. Figures 
|5]^a) and [5j^b) show the hysteretic spin-dependent voltages at some of the Vd’s among the 
multiple measurements. As shown in Figs|5](c) and [5j^d) , the magnitude of the resulting 
voltage change AV is maximal near the charge neutrality point (V G =13 V for SI and 
V G =10 V for S2). 

We do not observe any sign change in AU throughout the V G range, suggesting that with 
a fixed current bias, the carrier type change from n to p does not hip the direction of the 
measured spin polarization. This observation is fully consistent with dominant contributions 
to the spin polarization originating from the helical Dirac surface states since we use a dc 
current for the measurement. For a given crystal momentum, the spin polarization of holes is 
opposite to that of electrons. However, with a hxed DC current, the propagation direction of 
holes (below the Dirac point) is the reverse of the propagation of electrons (above the Dirac 
point). Therefore, the opposite spin polarization and the opposite propagation direction 
between electrons and holes make the sign of the detected spin-dependent voltages remain 
the same above and below the Dirac point. 12 We also note that our observations do not reveal 
any significant contributions from co-existing Rashba states which would have resulted in 
the opposite behavior due to the spin polarization being opposite to that of topological 
surface states 17 . For the latter, we assume that the Rashba bands minimum is close to the 
bottom of the conduction band and the chemical potential is tuned above the Rashba bands 
minimum. It is also possible that in the measurement range the chemical potential does not 
lie on any Rashba spin-split, bands. 
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F. Spin polarization calculation 


Equation (2) can be re-written as AV = IRbPtiPfm, where Pri is the spin polarization 
of the topological insulator and Pfm is the magnetic polarization of the ferromagnetic de¬ 
tector. The ballistic conductance <Jb representing the unit conductance times the number 
of propagating modes is given by: as = Rb~ 1 ~ where kp and w are the Fermi wave 

vector and the width of the topological insulator channel, respectively. For 3D bulk con¬ 
duction, kp = 37T 2 n where n is the 3D carrier density from all the bulk modes. However, 

with the chemical potential in the bulk bands, surface conduction may also co-exist along 
with the bulk conduction, with an associated kp = sJ\.Tvn s where n s is the 2D carrier density 
of a surface conduction channel. 32 Thus, a proper accounting of the distribution of current 
between bulk and surface conduction channels is actually needed to interpret experimental 
data when the chemical potential lies in the bulk bands. 

Our studies of ambipolar transport and weak anti-localization at different gate voltages 
suggest that, near the charge neutrality point, the bulk is depleted and a current flows 
only through the surface states. In this case, we can calculate kp using the lowest carrier 
density near the charge neutrality point (7.72 x 10 12 cm ~ 2 for Dev3 and 4.13 x 10 12 cm -2 
for Dev4); this corresponds to the residual carrier density from the puddles of electrons and 
holes created by disorder 32 33 . For the surface state, we determine kp = \J Anris = 0.070 A -1 
for Dev3 and 0.051 A -1 for Dev4. Using P FM fh 0.51 — 0.63 for a 20 nm thick Py layei®ES and 
the maximal AV near the charge neutrality point, we estimate the spin polarization of the 
topological insulator surface state Ppi = 0.42±0.15 for Dev3 and 0.78±0.26 for Dev4. These 
values are consistent within the error bars and comparable to theoretically calculated values 
of 0.5 36 W 2 / 7 r 17 . We note that the estimated numbers are lower bounds since the effective 
magnetic polarization P F m could have lower values for non-ideal A^CU/Py junctions. 


IV. CONCLUSIONS 

In conclusion, we have carried out electrical measurements of the spin polarization of 
the surface currents in electrically gated (Bi,Sb) 2 Te 3 devices as a function of the chemical 
potential. We estimate the spin polarization of the topological surface state to be 0.42 ±0.15 
and 0.78 ± 0.26 from two devices. Larger spin-dependent voltages were detected as the 










chemical potential approaches the Dirac point. Our results suggest that to achieve the 
maximum spin-dependent voltages from 3D topological insulators it is important to avoid 
the effects from the bulk bands, such as the parallel conduction of the unpolarized bulk 
carriers or the opposite spin polarization by Rashba spin-orbit interaction. Also, the results 
point toward a new strategy for electrical control of the magnitude of spin-dependent voltages 
from the current-induced spin polarization for potential topological spintronics, using gate- 
tunable topological insulator films. 
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TABLE I. List of devices with different MTJ areas, MTJ shapes, (Bi,Sb) 2 Te 3 samples, thicknesses 
of (Bi,Sb) 2 Te 3 film, and thicknesses of AI 2 O 3 tunnel barrier. 


Device MTJ Area 

(, fim 2 ) 

MTJ 

Shape 

Sample (Bi,Sb) 2 Te 3 

t (nm) 

AI 2 O 3 

t (nm) 

Devi 

50 2 7t 

circle 

SI 

7 

2.4 

Dev2 

100 2 tt 

circle 

SI 

7 

2.4 

Dev3 

50x50 

square 

SI 

7 

2.4 

Dev4 

10x80 

rectangle 

S2 

8 

2.1 

Dev5 

40x80 

rectangle 

S3 

6 

2.7 
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FIG. 1. (Color online) (a) Device schematic of an A^Os/Py MTJ on a cross-shaped topological 
insulator channel. Contact leads are labeled as 1 through 5. (b) A cross-section of an MTJ is 
illustrated, (c) An optical microscope image of a device with a circular-shaped MTJ of 50 /um 
radius (Devi). Si 3 N 4 layer is false-colored in green, (d) A cross-sectional HAADF-STEM image 
of the A^Oa/Py MTJ on (Bi,Sb) 2 Te 3 (Dev5), viewed along the [1 1 0] direction of the SrTiCH 
substrate, (e) Atomic-resolution HAADF-STEM image of the (Bi,Sb) 2 Te 3 structure, also viewed 
along the [1 1 0] direction of the substrate. STEM images were FFT low-pass filtered to 1.5 A(0.67 
A -1 in reciprocal space) resolution to remove instrumental noise. (Bi,Sb) atoms are indicated in 
blue, Te atoms are indicated in yellow, and 1 quintuple layer (QL) is labeled. 


13 

























JJ 


10 


S3 


FIG. 2. (Color online) (a) Gate voltage dependence of the two-terminal channel resistance (red 
circles) at zero field and Hall resistance (green triangles) at H = 10 kOe from Devi (SI). Inset: 
measurement setup, (b) Similar measurements as (a) for Dev4 (S2). Inset: schematic of the 
band structure of (Bi,Sb) 2 Te 3 in presence of Vq- (c) Inverse of the Hall resistance for Devi (green 
triangles) and for Dev4 (blue circles) in units of 2D carrier density (cm -2 ). Inset: zoom-in near the 
charge neutrality point (Vq = 13 V for Devi and Vq = 10 V for Dev4). (d) Variation of prefactor 
a with Vg, obtained by fitting magneto-resistance (inset) to Eq. (1) for Devi (green triangles) 
and for Dev4 (blue circles). Error bars represent an uncertainty with 95 percent confidence in data 
fitting. Inset: magneto-resistance for Dev4 at Vq = —140 V (pink), —20 V (blue), 20 V (cyan), 60 
V (orange), 140 V (red). All data from Devi (Dev4) were taken at 2.1 K (1.8 K). 
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FIG. 3. (Color online) (a) Temperature dependence of the RA product shows an insulating behav¬ 
ior. Bottom inset illustrates a schematic of three-terminal measurement setup. Top inset shows 
temperature dependence of the channel resistance of Devi (SI) and Dev4 (S 2 ). (b )Vdc depen¬ 
dence of RA product from Devi at 1.8 K with different gate voltages. Middle inset illustrates the 
tunneling across the AI 2 O 3 tunnel barrier between Voc'-tuned ferromagnetic metal and topological 
insulator top surface. 
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FIG. 4. (Color online) Detected spin-dependent voltages with positive DC bias current (a) from 
Dev4 [3 /iV (red) and 1 fi\ (green)] and (b) from Dev2 [30 /rV (red), 10 fAT (gray), and 3 /rV 
(green)]. Reversed spin-dependent voltages with negative bias currents (c) from Dev4 [—4 (blue) 
and —1 /rV (cyan)] and (d) from Dev2 [—30 /rV (blue)]. The schematics depict spin polarization (red 
arrow) induced by electron current (black arrow) on a surface of a topological insulator channel and 
the magnetization of the ferromagnetic layer (white arrow). The voltage change AV as a function 
of bias current (e) from Dev4 and (f) from Dev2 shows a linear relationship. The black line is a 
linear fit to data shown as red circles with error bars that represent the standard deviations over 
multiple measurements. Inset in (e) illustrates the measurement setup. All data were taken at 1.8 
K with Vg= 0 V. 
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FIG. 5. (Color online) Detected hysteretic spin-dependent voltages with up sweep (red) and down 
sweep (black) of in-plane magnetic field with different gate voltages (a) from Dev3 and (b) from 
Dev4. Observed AV with respect to Vq is plotted in (c) for Dev3 and in (d) for Dev4. The error 
bars represent the standard deviations over multiple measurements. All data were taken with fixed 
currents of 30 fiA for Dev3 and 2 jj ,A for Dev4 at 1.8 K. 
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